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Using recently published, high-precision 7r"'"7r~ cross section data by the BABAR experiment from 
the analysis of e+e" events with high-energy photon radiation in the initial state, we reevaluate the 
lowest order hadronic contribution (ajl^'^''^'^) to the anomalous magnetic moment of the muon. 
We employ newly developed software featuring improved data interpolation and averaging, more 
accurate error propagation and systematic validation. With the new data, the discrepancy between 
the e^e~ and r-based results for the dominant two-pion mode reduces from previously 2Aa to 1. Sa- 
in the dispersion integral, though significant local discrepancies in the spectra persist. We obtain for 
the e"'"e~-based evaluation aji^'*''^° = (695.5 ± 4.1) • 10"^^°, where the error accounts for all sources. 
The full Standard Model prediction of differs from the experimental value by 3.2a. 



I. INTRODUCTION 

The Standard Model (SM) prediction of the anomalous 
magnetic moment of the muon, a^, is limited in preci- 
sion by contributions from hadronic vacuum polarisation 
(VP) loops. These contributions can be conveniently sep- 
arated into a dominant lowest order {a^^'^'^) and higher 
order (ajj'*'^'^^) parts. The lowest order term can be cal- 
culated with a combination of experimental cross section 
data involving e+e~ annihilation to hadrons, and per- 
turbative QCD. These are used to evaluate an energy- 
squared dispersion integral, ranging from the tt'^j thresh- 
old to infinity. The integration kernel strongly empha- 
sises the low-energy part of the spectrum. About 73% of 
the lowest order hadronic contribution is provided by the 
7r+7r^(7) final state. "'^ More importantly, 62% of its to- 
tal quadratic error stems from the tt+tt" mode, stressing 
the need for ever more precise experimental data in this 
channel to confirm or not the observed deviation of 3.8cr 
between SM prediction and experiment [1]. 

A former lack of precision e"'"e~ -annihilation data in- 
spired the search for an alternative. It was found [2] in 
form of T ^ Vr + TT~n^ , 27r~7r+7r'', tt^Stt^ spectral func- 
tions [3-6], transferred from the charged to the neutral 
state using isospin symmetry. During the last decade, 
new measurements of the tt^tt" spectral function in e'^e~ 
annihilation with percent accuracy became available [7- 
10] , superseding or complementing older and less precise 
data. With the increasing precision, which today is on a 
level with the r data in that channel, systematic discrep- 
ancies in shape and normalisation of the spectral func- 
tions were observed between the two systems [11, 12]. It 
was found that, when computing the hadronic VP contri- 
bution to the muon magnetic anomaly using the r instead 



^ Throughout this paper, final state photon radiation is implied 
for the 7r+7r~ final state. 



of the e+e~ data for the 27r and An channels, the observed 
deviation with the experimental value [13] would reduce 
to less than la [14]. 

The discrepancy with the e+e" data decreased af- 
ter the inclusion of new r data from the Belle exper- 
iment [15], published e+e~ data from CMD2 [8] and 
KLOE [10] (superseding earlier data [16]), and a reevalu- 
ation of isospin-brcaking corrections affecting the r-based 
evaluation [1].-^ In terms of a^'^'^'^ , the difference be- 
tween the r and e+e~-based evaluations in the dominant 
7r+7r~ channel was found to be 11.7 ± 3. See i 3.5^+18 [1] 
(if not otherwise stated, this and all following num- 
bers are given in units of 10~^°), where KLOE exhibits 
the strongest discrepancy with the r data (without the 
KLOE data the discrepancy reduces from 2.4cr to 1.9cr). 
Another quantity for comparison, which is more sensitive 
to the higher-energy tt+tt" spectrum, is the t~ ■k~-k^v 
branching fraction showing a difference between measure- 
ment and e+e" prediction of (0.64±0.10T-±0.28ee)% [1].^ 

Recently, the BABAR Collaboration has published [17] 
a 7r+7r~(7) spectral function measurement based on half 
a million selected e+e" — )• 7r"'"7r~7(7) events, where the 
hard photon is dominantly radiated in the initial state 
(ISR) . It benefits from a large cancellation of systematic 
effects in the ratio 7r+7r~7(7) to /i+/i~7(7) employed for 
the measurement. In this letter, wc present a reevalua- 
tion of the lowest order hadronic contribution to in- 
cluding the new BABAR data. We deploy the new soft- 
ware package HVPTools [18], featuring a more accurate 
data interpolation, averaging and integration method, 



^ The total size of the isospin-breaking correction to a^^'^'^'~' has 
been estimated to (—16.1 ± 1.9) ■ 10"-'^'', which is dominated by 
the short-distance contribution of (-12.2 ± 0.2) ■ 10"^" [1]. 

^ A total isospin-breaking correction of (-1-0.69 ± 0.22)% has been 
added to the e+e" prediction of the T~ — > tt^tt"!/ branching 
fraction [1]. 
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better systematic tests, and an improved statistical anal- 
ysis based on the generation of large samples of pseudo 
experiments. 

II. e+e-^ir+7T- CROSS SECTION DATA 

The dispersion integral for the lowest order hadronic 
contribution reads 

1 f°° 

^had.LO ^ / dsK(s)ae+e-^hadrons(s) , (1) 

where K(s) ~ s^^ [19]. The contribution from the light 
u, d, s quark states is evaluated using exclusive experi- 
mental cross section data up to an energy of 1.8 GeV, 
where resonances dominate, and perturbative QCD to 
predict the quark continuum beyond that energy. In 
this work we only reevaluate the contributions from the 
e+e~ — )• TT+TT" and 7r+7r~27r° channels. For all the others 
wc refer to Refs. [11, 12, 14]. 

A large number of e+e~ — >■ tt+tt" cross section mea- 
surements are available. Older measurements stem from 
OLYA [20, 21], TOF [22], CMD [20], DM1 [23] and 
DM2 [24].'* They are affected by an incomplete or un- 
documented application of radiative corrections. Equa- 
tion (1) and the treatment of higher order hadronic con- 
tributions require initial state radiation as well as lep- 
tonic and hadronic VP contributions to be subtracted 
from the measured cross section data, while final state 
radiation should be included. Because of lack of docu- 
mentation, the latter contribution of approximately 0.9% 
in the tt+tt" channel has been added to the data, ac- 
companied by a 100% systematic error [11]. Initial state 
radiation and leptonic VP effects are corrected by all ex- 
periments, however hadronic VP effects are not. They 
are strongly energy dependent, and in average amount 
to approximately 0.6%. We apply this correction accom- 
panied by a 50% systematic error [11]. These FSR and 
hadronic VP systematic errors are treated as fully corre- 
lated between all measurements of one experiment, and 
also among different experiments. 

More recent precision data, where all required radia- 
tive corrections have been applied by the experiments, 
stem from the CMD2 [8] and SND [9] experiments at the 
VEPP-2M collider (Novosibirsk, Russia). They achieve 
comparable statistical errors, and energy-dependent sys- 
tematic uncertainties down to 0.8% and 1.3%, respec- 
tively. 

These measurements have been complemented by re- 
sults from KLOE [10] at DA$NE (Frascati, Italy) run- 
ning at the <j) resonance centre-of-mass energy. KLOE 
applied for the first time a hard-photon ISR technique to 
precisely determine the tt+tt" cross section between 0.592 



and 0.975 GeV. The cross section data are obtained from 
a binned distribution, corrected for detector resolution 
and acceptance effects. The analysed data sample cor- 
responds to 240 pb~* integrated luminosity providing a 
0.2% relative statistical error on the tt+tt" contribution 
to a^'^'^'^. KLOE does not normalise the 7r+7r~7 cross 
section to e+e~ ;U+/i~7 so that the ISR radiator func- 
tion must be taken from Monte Carlo simulation (c/. [26] 
and references therein). The systematic error assigned to 
this correction varies between 0.5% and 0.9% (closer to 
the (j) peak) . The total assigned systematic error lies be- 
tween 0.8% and 1.2%. 

In a recent publication [17] the BABAR Collabora- 
tion reported measurements of the processes e+e~ — )• 
7r+7r~ (7), /i^/i" (7) using the ISR method at 10.6 GeV 
centre-of-mass energy. The detection of the hard ISR 
photon allows BABAR to cover a large energy range 
from threshold up to 3 GeV for the two processes. The 
7r'^7r^(7) cross section is obtained from the 7r+7r~7(7) 
to /U+/(x~7(7) ratio, so that the ISR radiation function 
cancels, as well as additional ISR radiative effects. Since 
FSR photons are also detected, there is no additional 
uiiccatainty from radiative corrections at NLO level. Ex- 
perimental systematic uncertainties are kept to 0.5% in 
the p peak region (0.6-0.9 GeV), increasing to 1% out- 
side. 



III. COMBINING CROSS SECTION DATA 

The requirements for averaging and integrating cross 
section data are: (i) properly propagate all the uncertain- 
ties in the data to the final integral error, {ii) minimise 
biases, i.e., reproduce the true integral as closely as pos- 
sible in average and measure the remaining systematic 
error, and {Hi) optimise the integral error after averag- 
ing while respecting the two previous requirements. The 
first item practically requires the use of pseudo-Monte 
Carlo (MC) simulation, which needs to be a faithful rep- 
resentation of the measurement ensemble and to contain 
the full data treatment chain (interpolation, averaging, 
integration) . The second item requires a flexible data in- 
terpolation method (the trapezoidal rule is not sufficient 
as shown below) and a realistic truth model used to test 
the accuracy of the integral computation with pseudo- 
MC experiments. Finally, the third item requires optimal 
data averaging taking into account all known correlations 
to minimise the spread in the integral measured from the 
pseudo-MC sample. 

The combination and integration of the e+e" — >■ tt+tt" 
cross section data is performed using the newly devel- 
oped software package HVPTools [18].^ It transforms 



" We do not use the data from NA7 [25]. 



^ HVPTools is written in object-oriented C++ and relies on ROOT 
functionality [27]. The cross section database is provided in XML 
format and can be made available to users - please contact the 
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the bare cross section data and associated statistical and 
systematic covariance matrices into fine-grained energy 
bins, taking into account to our best knowledge the cor- 
relations within each experiment as well as between the 
experiments (such as uncertainties in radiative correc- 
tions). The covariance matrices are obtained by assum- 
ing common systematic error sources to be fully corre- 
lated. To these matrices are added statistical covariances, 
present for example in binned measurements as provided 
by KLOE, BABAR or the t data, which are subject to 
bin-to-bin migration that has been unfolded by the ex- 
periments, thus introducing correlations. 

The interpolation between adjacent measurements of 
a given experiment uses second order polynomials. This 
is an improvement with respect to the previously applied 
trapezoidal rule, corresponding to a linear interpolation, 
which leads to systematic biases in the integral (see be- 
low, and also the discussion in Sec. 8.2 and Fig. 12 of 
Ref. [11]). In the case of binned data, the interpolation 
function within a bin is renormalised to keep the integral 
in that bin invariant after the interpolation. This may 
lead to small discontinuities in the interpolation function 
across bin boundaries. The final interpolation function 
per experiment within its applicable energy domain is 
discretised into small (1 MeV) bins for the purpose of av- 
eraging and numerical integration. 

The averaging of the interpolated measurements from 
different experiments contributing to a given energy bin 
is the most delicate step in the analysis chain. Corre- 
lations between measurements and experiments must be 
taken into account. Moreover, the experiments have dif- 
ferent measurement densities or bin widths within a given 
energy interval and one must avoid that missing informa- 
tion in case of a lowc;r measurement density is substituted 
by extrapolated information from the polynomial inter- 
polation. To derive proper averaging weights given to 
each experiment, wider averaging regions^ are defined to 
ensure that all locally available experiments contribute 
to the averaging region, and that in case of binned mea- 
surements (KLOE, BABAR, t data) at least one full bin 
is contained in it. The averaging regions are used to com- 
pute weights for each experiment, which are applied in 
the bin- wise average of the original finely binned interpo- 



authors. The systematic errors are introduced component by 
component as an algebraic function of mass or as a numerical 
value for each data point (or bin). Systematic errors belonging 
to the same identifier (name) are taken to be fully correlated 
throughout all measurements affected. So far, HVPTools has 
been only employed for the numerical evaluation of the most 
important tt+tt" (and 7r+7r~27r'') parts of the dispersion inte- 
gral (1). 

For example, when averaging two binned measurements with un- 
equal bin widths, a useful averaging region would be defined by 
the experiment with the larger bin width, and the bins of the 
other experiments would be statistically merged before comput- 
ing the averaging weights. 



lation functions.^ If the value of a bin- wise average* 
exceeds the number of degrees of freedom (ridof), the error 
in this averaged bin is rescaled by \/x^ /n^oi to account 
for inconsistencies (c/. Fig. 1). Such inconsistencies fre- 
quently occur because most experiments are dominated 
by systematic uncertainties, which are difficult to esti- 
mate. 

The consistent propagation of all errors into the eval- 
uation of a}^'^'^'^ is ensured by generating large samples 
of pseudo experiments, representing the full list of avail- 
able measurements and taking into account all known 
correlations. For each generated set of pseudo measure- 
ments, the identical interpolation and averaging treat- 
ment leading to the computation of Eq. (1) as for real 
data is performed, hence resulting in a probability den- 
sity distribution for a^'"^'^'^ (tt+tt"), the mean and RMS 
of which define the la allowed interval (and which - by 
construction - has a proper pull behaviour). The pro- 
cedure yielding the weights of the experiments can be 
optimised with respect to the resulting error on a}^'^'^'^ . 

We have tested the fidelity of the full analysis chain 

(polynomial interpolation, averaging, integration) by 
using as truth representation a Gounaris-Sakurai [28] 
vector-meson resonance model faithfully describing the 
tt+tt" data. The central values for each of the avail- 



The averaging weights for each experiment are computed as fol- 
lows: 

1. pseudo-MC generation fluctuates the data points (or bins) 
along the original measurements taking into account all 
known correlations; the polynomial interpolation is redone 
for each generated pseudo MC; 

2. the averaging regions are filled for each experiment and each 
pseudo-MC generation and interpolated with second order 

polynomials; 

3. small (1 MeV) bins are filled for each experiment, in the 
energy intervals covered by that experiment, using the in- 
terpolation of the averaging regions; 

4. in each small bin a correlation matrix between the exper- 
iments is computed from which the averaging weights are 

obtained. 

The bin-wise average between experiments is computed as fol- 
lows: 

1. pseudo-MC generation fluctuates the data points (or bins) 
along the original measurements taking into ax;count all 
known correlations; the polynomial interpolation is redone 

for each generated MC; 

2. for each generated pseudo-MC, small (1 MeV) bins are filled 
for each experiment, in the energy intervals covered by that 
experiment, using the polynomial interpolation; 

3. the average and its error arc computed in each small bin 
using the weights previously obtained; 

4. the covariance matrix among the experiments is computed 

in each small bin; 

5. rescaling corrections are computed and applied for each 
bin. 
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FIG. 1: Rescahng factor accounting for inconsistencies among 
experiments versus ^/s (see text). The peak around O.4 GeV 
is introduced by a discrepancy between CMD2 and TOF 
measurements versus BABAR. The peaks around 0.65 and 
0. 74 GeV are introduced by outlier from CMD. The sharp peak 
at 0. 78 Ge V is due to local discrepancies along the p-oj in- 
terference. The bump between 0.85 and 0.95 GeV is due to 
a discrepancy between KLOE and BABAR. Finally, between 
1.4.5 and 1.65 GeV measurements from MEA and DM2 sig- 
nificantly exceed the BABAR cross sections. 



able measurements are shifted to agree with the Breit- 
Wigner model, leaving their statistical and systematic 
errors unchanged. The so created set of measurements 
is then analysed akin to the original data sets. The dif- 
ference between true and estimated a}},'^'^'^ values is a 
measure for the systematic uncertainty due to the data 
treatment. Wc find negligible bias below 0.1 (remember 
the 10~^° unit), increasing to 0.5 (1.2 without the high- 
density BABAR data) when using the trapezoidal rule 
for interpolation instead of second order polynomials. 

The individual e~^e~ — >■ tt+tt" cross section measure- 
ments (dots) and their average (shaded/green band) are 
plotted in Fig. 2. The error bars contain statistical and 
systematic errors. For better comparison wc also plot in 
Fig. 3 the relative differences between BABAR, KLOE, 
CMD2, SND, and the average. Fair agreement is ob- 
served, though with a tendency to larger (smaller) cross 
sections above -0.8 GeV for BABAR (KLOE). These in- 
consistencies (among others) lead to the error rescaling 
shown versus ^y.s in Fig. 1. 

The left hand plot of Fig. 4 shows the weights ver- 
sus ^/s the different experiments carry in the average. 
BABAR and KLOE dominate over the entire energy 
range. Owing to the sharp radiator function, the avail- 
able statistics for KLOE increases towards the (p mass, 
hence outperforming BABAR above ~0.8GeV. For ex- 
ample, at 0.9 GeV KLOE data have statistical errors of 
0.5%, which is twice smaller than for BABAR (renor- 
malising BABAR to the 2.75 times larger KLOE bins 
at that energy). Conversely, at 0.6 GeV the compar- 
ison reads 1.2% (KLOE) versus 0.5% (BABAR, again 
given in KLOE bins which are about 4.2 times larger 



than BABAR at that energy). The experiments labelled 
"other exp" in the figure correspond to older data with 
incomplete radiative corrections. Their weights are small 
throughout the entire energy domain. 

Figure 4 (right) shows versus ^/s the combined e+e^ — )■ 
TT^TT" cross section multiplied by the kernel function 
K{s) occurring in the dispersion integral (1). The ker- 
nel strongly emphasises the low-energy spectrum. The 
dashed (red) curve belonging to the right axis in the 
plot gives the corresponding error contribution (diagonal 
errors only, statistical and systematic errors have been 
added in quadrature). The peaks arc introduced by the 
error rescaling and indicate inconsistencies between the 
measurements. The uncertainty in the integral is domi- 
nated by the measurements below 0.8 GeV. 



IV. RESULTS 

A compilation of results for aj^'^'^'^^ [tttt] for the var- 
ious sets of experiments and energy regions is given in 
Table L The comparison with our previous result [1], 
^had,LOj^^j ^ 503 5 _|_ 3.5^.^^^ shows that the inclusion 

of the new BABAR data significantly increases the cen- 
tral value of the integral, without however providing a 
large error reduction. This is due to the incompati- 
bility between mainly BABAR and KLOE, causing an 
increase of the combined error. In the energy interval 
between 0.63 and 0.958 GeV, the discrepancy between 
evaluations from KLOE and BABAR 
amounts to 2.0tT. BABAR is the only experiment cov- 
ering the entire energy region between 2771,^ and 1.8 GeV. 
Using only the BABAR data to evaluate 

^had,LO[^^j one 

finds [17] 514.1 ± 2.2,tat ± 3.1,y,f 

Also given in Table I is the combined r-based result 
from Ref. [1]. The difference between the t and e~^e~- 
based evaluations of ajj'*'^'^*-'[7r7r] now reads 6.8±3.5r+iB± 
2.9ee, thus reducing to 1.5a compared to 2.4tT without 
BABAR [1] (the BABAR-only result is in excellent agree- 
ment with the r data).® A comparison between the com- 
bined e^e~ and t two-pion cross sections relative to the 
e+e" result is shown in Fig. 5. Significant local discrep- 
ancies arise in particular above the p peak. 

We also reevaluate the e+e" — > 7r"'"7r~27r*' contribu- 
tion to a^^'^'^'^. The CMD2 data used previously [31] 
have been superseded by modified or more recent, but 
yet unpublished data [32], recovering agreement with 
the published SND cross sections [33]. Since the new 
data are unavailable, we discard the obsolete CMD2 data 
from the 7r+7r~27r° average, finding a^'^'^'^'^ [7:1:2-1:°] = 
17.6±0.4gtat±l-7syst (compared to 17.0±0.4stat±l-6syst 



^ Combining the e+e^ and r-bascd evaluations would lead to an 
error rescaling by a factor of 1.5 to account for the inconsistency 
in the integrated data. This would approximately cancel the 
expected precision gain from the combination. 
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FIG. 2: Cross section for e^e~ — > tt+tt^ annihilation measured by the different experiments for the entire energy range (top), 
and zoomed energy intervals (all other plots). The errors bars contain both statistical and systematic errors, added in quadrature. 
The shaded (green) band represents the average of all the measurements obtained by HVPTools, which is used for the numerical 
integration following the procedure discussed in Sec. III. 



6 




0.9 1 
^Js [GeV] 



0.2 ft 

■ 

g> 0.15 



0.1 



> 
< 



^ 0.05 

X 

^ -0 
o 

1 -0.05 
a> 

(0 

« -0.1 
o 

O -0.15 



■o-^o'^r 



I ' ' ' ' I — ' — ^ 

Average 
' KLOE 



0.6 



0.7 



0.8 



0.9 1 
\Js [GeV] 



'0.5 




0.6 



0.7 



0.8 



Average 
CMD2 06 
CMD2 03 

0.9 1 
nIs [GeV] 



0.2 
g> 0.15 



0.1 



> 
< 



2 0.05 

X 

o 

Z -0.05 
a> 

(A 

« -0.1 
o 

O -0.15 



■°-2oV 




0.6 



0.7 



0.8 



Average 
SND 



0.9 1 
vjs [GeV] 



FIG. 3: Relative cross section comparison between individual experiments (symbols) and the HVPTools average (shaded band) 
computed from all measurements considered. Shown are BABAR (top left), KLOE (top right), CMD2 (bottom left) and SND 
(bottom right). 
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FIG. 4: Left: relative averaging weights per experiment versus ^/s. Right: contribution to the dispersion integral (1) for the 
combined e^ e~ data obtained by multiplying the -k'^-k~ cross section by the kernel function K{s) (solid line). The dashed (red) 
curve belonging to the right axis shows the corresponding error contribution, where statistical and systematic errors have been 
added in guadrature. Note that the information conveyed by this curve is incomplete because only diagonal errors are shown, 
disregarding correlations between the cross section measurements which have significant influence on the integral error. 
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TABLE I: Evaluated a|i'"*''^ [tttt] contributions from the e"'"e~ data for different energy intervals and experiments. Where two 
errors are given, the first is statistical and the second systematic. Also given is the r-based result from Ref. [1] combining all 
available t data. The combined error has been rescaled to account for the inconsistency between the two evaluations. 



Energy range (GeV) 


Experiment 




2m^± — 0.3 


Combined e^e~ (fit) 


0.55 ± 0.01 


n on A ^; o 
U.oU — U.DO 


Combined e'''e~ 


ioz.D ± U.8 ± i.U (i.Otot j 


u.Do — u.yoo 




1 O _i_ O A _J_ O 1 / Q O ^ 

ODi.O ± Z.4 ± i.i 




CMD2 06 


360.2 ± 1.8 ± 2.8 (3.3tot) 




SND 06 


360.7 ± 1.4 ±4.7 (4.9tot) 




KLOE 08 


356.8 ±0.4 ±3.1 (3.1tot) 




BABAR 09 


365.2 ± 1.9 ± 1.9 (2.7tot) 




Combined e^e~ 


360.8 ± 0.9 ± 1.8 (2.0tot) 


0.958- 1.8 


Combined e^e~ 


14.4 ±0.1 ±0.1 (0.2tot) 


Total 


Combined e^e~ 


508.4 ± 1.3 ±2.6 (2.9tot) 


Total 


Combined r [1] 


515.2 ± 2.0exp ± 2.28 ± l-9iB (3.5tot) 




FIG. 5: Relative comparison between the combined r (dark 
shaded) and e'^e~ spectral functions (light shaded), nor- 
malised to the e^e~ result. The apparently oscillating struc- 
ture around 0.5 GeV is due to two Belle measurements fluc- 
tuating to large cross section values. Clearly visible is the 
interference due to p-(j) mixing around 1 GeV, which is not 
included in the is o spin-breaking corrections applied to the r 
data. It is also visible in the upper, and lower nght hand plots 
of Fig. 2. The deviation between 0.8 and 0.95 GeV is due to 
the discrepancy between t and KLOE data, which dominate 
in this region (cf. Fig. 4 left)- Comparing the r data with 
the combined e^e~ data instead of a fit to a single experiment 
CMD-2 limited to 1 GeV as it was done for Fig. 4 in Ref- [29] 
and Fig. 28 m Ref. [30], we observe a reduced discrepancy, 
in particular between l.OGeV and 1.4GeV. We therefore dis- 
agree with the conclusion reached m these references, where 
the difference goes up to a factor 4, and ts even in the oppo- 
site direction with respect to the one we observe. 



FIG. 6: Cross section measurements e e~ — > vr 7r~27r used 
in the calculation of a)]^'^'^'~'[iT-K2-K'^]. The shaded band depicts 
the HVPTools interpolated average within la errors. The in- 
dividual measurements are referenced in [11]. 

when including the obsolete CMD2 data). The corre- 
sponding cross section measurements and HVPTools av- 
erage are shown in Fig. 6. 

Adding to the e+e"-based o^^^-^'^I'ktt] and 
^had,LO |-.^^27r°] results the remaining exclusive multi- 
hadron channels as well as perturbative QCD [14], we 
find for the complete lowest order hadronic term 

^had,LO[g+g-] ^ 5 ^ 4_o^^p ± o.7qcd (4.1tot) ■ 

It is noticeable that the error from the tt+tt" channel now 
equals the one from all other contributions to a^'^'^'^ . 

Adding further the contributions from higher order 
hadronic loops, —9.79 ± O.OSoxp ± O.OSrad [34], hadronic 
light-by-light scattering (LBLS), 10.5 ± 2.6 [35], as well 
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-285 ±51 1- 

JN 09 (e*el 

-299 ±65 I 

Davier et al. 09 (x) 
-157 ±52 

Davier et al. 09 (e*e") 
-312±51 I 

This work (eV w/ BABAR) 
-255 ±49 

BNL-E821 (WA) 

0±63 
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FIG. 7: Compilation of recent results for , subtracted 
by the central value of the experimental average [40]. The 
shaded vertical band indicates the experimental error. The 
SM predictions are taken from: HMNT 07 [34], JN 09 [30], 
Davier et al. 09 [1] (r-based and e'^e" including KLOE), and 
the e"'"e~-based value from this work. 



as QED, 11 658 471.809± 0.015 [36] (see also [37] and ref- 
erences therein), and electroweak effects, 15.4 ± O.lhad i 
0.2Higgs [38, 39], we obtain the SM prediction (still in 
10-1° units) 

a™[e+e"] = 11 659 183.4 ± 4.1 ± 2.6 ± 0.2 (4.9tot) , 

where the errors have been split into lowest and higher or- 
der hadronic, and other contributions, respectively. The 
a^[e^e~] value deviates from the experimental average, 
a^'^P = 11 659 208.9±5.4±3.3 [37, 40], by 25.5±8.0 (3.2cr). 

A compilation of recent SM predictions for com- 
pared with the experimental result is given in Fig. 7. 
The BABAR results are not yet contained in evaluations 
preceding the present one. The result by HMNT [34] con- 
tains older KLOE data [16] , which have been superseded 
by more recent results [10], leading to a slightly larger 
value for a]','^'^'^^ . 



V. CONCLUSIONS 

We have reevaluated the lowest order hadronic contri- 
bution to the niuon magnetic anomaly in the dominant 
7r+7r~ channel, using new precision data published by 
the BABAR Collaboration. After combination with the 
other e+e^ data a 1.5cr difference with the r data re- 
mains for the dominant 7r~'"7r~ contribution. For the full 
e"'"e~ -based Standard Model prediction, including also 
a reevaluated 7r+7r~27r° contribution, we find a devia- 
tion of 3.2(T from experiment (reduced from 3.7(T with- 
out BABAR). The deviation reduces to 2.9cr when ex- 
cluding KLOE data, and further decreases to 2.4(t when 
using only the BABAR data in the tt+tt" channel. As 
a reminder, the r-based result deviates by 1.9(7 from the 
Standard Model. 

The present situation for the evaluation of aj^'*'i''"'-'[7r7r] 
is improved compared to that of recent years, as more in- 
put data from quite different experimental facilities and 
conditions have become available: e~^e~ energy scan, 
e+e^ ISR from low and high energies, t decays. Our 
attitude has been to combine all the data and include in 
the uncertainty the effects from differences in the spec- 
tra. At the moment the ideal accuracy cannot be reached 
as a consequence of the existing discrepancies due to un- 
corrected or unaccounted systematic effects in the data. 
A critical look must be given to the different analyses 
in order to identify their weak points and to improve on 
them or to assign larger systematic errors. 

It is thereby not sufficient to concentrate on improv- 
ing the 7r"'"7r~ channel alone. Problems also persist in the 
7r+7r^27r° mode, where the t and e^e~ -based evaluations 
differ by (3.8±2.2) • 10"^°, but also the e'^e~ data among 
themselves exhibit discrepancies. Fortunately, new pre- 
cision data from BABAR should soon help to clarify the 
situation in that channel. 
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